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A review of discrete quantum walk with two particle is given. The use of different states encoun- 
tered in identical particle, and the idea of entanglement and superposition is explored to explored 
the interesting dynamics of two particle quantum walk. Boundary conditions can specify certain 
dynamics and so a survey of periodic boundary condition (circle) is presented. A simulation for a 
Hadamard walk for different periods of a circle is considered and results are compared for various 
periods and for a case with absorbing boundaries. 



A. Introduction 

Quantum information is represented in the form of 
qubits which are similar to classical bits but have greater 
information carrying and processing capacities. Thus 
as a result superposition of states under quantum dy- 
namics, which progress by exploring multiple possible 
paths simultaneously with the amplitudes correspond- 
ing to different paths interfering. Some of the algo- 
rithms in quantum information has been shown to have 
a successful implementation using quantum walks mo- 
tivating why we study quantum walk in this piece of 
work. Quantum walks was introduced in 199cP to repre- 
sent a quantum version of the classical walk, and has be- 
come an active field of research in quantum information 
and computation^. The variance of quantum walk on a 
line grows quadratically with the number of steps (time) , 
compared to the linear growth for the classical random 
wallM A probabilistic result is obtained in quantum walk 
upon measurement. Quantum walk has been broadly 
studied in two different forms: continuous-time quantum 
wallP and discrete-time quantum walkP. A HadamarcP 
discrete-time quantum walk is considered in this work. 
A background tool is presented to capture the nature of 
dynamics. Results of two particle on a circle is presented 
and it comparison with a straight line is made. 



I. DYNAMICS OF QUANTUM WALK 

Quantum dynamics show characteristics that do not 
have a classical analogue but in many cases of studies 
just like their classical analogue, they are considered in 
parallel. This is a similar case of quantum walk which is 
widely studied in two forms: continuous-time quantum 
walk and discrete-time quantum walk. In continuous- 
time quantum walk, there is a direct definition of the 
walk on the position space^, whereas in the discrete-time 
quantum walk, there is a necessity to introduce a quan- 
tum coin operation to define the direction in which the 
particle has to move^. These two approaches look re- 
ally different but the results are often similar, but due to 
the coin degree of freedom, the discrete-time variant has 



been shown to be more powerful than the other in some 
contexiP. To match the performance of the discrete-time 
quantum walk, the coin degree of freedom can be intro- 
duced in the continuous-time quantum wallP. Discrete- 
time quantum walk is the point of concern for our study. 



A. Continuous-time quantum walk 

Similar to a technique usually used in physics, it is 
possible to define the continuous-time quantum walk, as 
the quantization of a classical walk. This can be achieved 
by introducing quantum amplitudes in place of classical 
probabilities . To do this it is easier to first define the 
continuous-time classical random walk. 

The continuous-time classical random walk takes place 
entirely in the position space. Consider a continuous- 
time classical random walk on the position space Hp 
spanned by a vertex set V of a graph G with edges set 
E, G = (V,E). An adjacency matrix is used for a step 
of the random walk A which transform the probability 
distribution over V, i.e., 



A 



j,k 



1 (j, k) £ E 
(j, k)$E 



(1) 



for every pair j, k € V. The matrix G is the generator 
matrix H given by 



d ]l 3 



H 







(j,k) e E 

otherwise 



(2) 



where dj is the degree of the vertex j and 7 is the tran- 
sition probability between neighboring nodes at any in- 
stant of time. 

Let Pj(t) be the probability of being at vertex j at 
time t then we can take the transition on graph G to be 
the solution of differential equation 

d ;P j (t) = -J2^,kPk(t). (3) 

kev 



dt 



with solutions given by 
P(t) 



HI 



P(0). 



(4) 
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By using quantization technique in which the probabil- 
ities Pj are replaced by quantum amplitudes a,j(t) — 
(j\ip(t)} where \j) is spanned by the orthogonal basis of 
the position Hilbcrt space T~i v and introducing a factor of 
i we obtain 



" dt 



kev 



We can see that ^ is the Schrodinger equation 
. d , 



' dV 



HI 



(5) 



(6) 



Since generator matrix is a Hermitian operator, the nor- 
malization is preserved during the dynamics. The so- 
lution of the differential equation can be written in the 
form 



\m) 



sHt 



W0)>. 



(7) 



Therefore, the continuous-time quantum walk is of the 
form of Schrodinger equation, a non-relativistic quantum 
evolution. 

In the implementation of continuous-time quantum 
walk on a line, the position Hilbert space % p can be writ- 
ten as a state span-fj^)}, where jeZ. The Hamiltonian 
H is defined such that, 



(8) 



and is made to evolve with time t by applying the trans- 
formation 

U(t) = exp(-iHt), (9) 
for time independent Hamiltonian H , and 



U(t) = exp(-i J Hdi), 



(10) 



for the case where the Hamiltonian is time dependent. 
The Hamiltonian H of the process acts as the generator 
matrix which will transform the probability amplitude at 
the rate of 7 to the neighboring sites, where 7 is time- 
independent constant. 



B. Discrete-time quantum walk 

We will first define the structure of the discrete-time 
classical random walk. The discrete-time classical ran- 
dom walk takes place on the position Hilbert space Hp 
with instruction from the coin operation. The coin could 
be fair or biased. A coin flip is used to define the direction 
in which the particle moves and a shift operation in suc- 
cession to moves the particle in position space depending 
on the outcome of the coin. A simple walk on a line, a 
fair coin with head and tail can define the movement to 
the left and right respectively. 



The discrete-time quantum walk also has a very similar 
structure to that of its classical counterpart. The relevant 
degrees of freedom are the particle's (walker) position i 
(with igZ) on the line, as well as its coin state. We can 
then define the elements of a discrete quantum walk as 
follows: 

Particle or walker: A quantum system living in 
a Hilbert space of infinite but countable dimension Hp 
which is spanned by {\i}} ■ The walker will be initialized 
at the origin without lost of generality. 

Coin: A quantum system living in a 2-dimensional 
Hilbert space T-L c with span span{ | t)> I 4-) } ■ This will 
be the quantum equivalent of randomly choosing which 
way the particle will move (like tossing a coin in the clas- 
sical case). The initial coin state will depends on the 
symmetry we want to imprint on the position probabil- 
ity distribution of the walker. 

The total Hilbert space is given by "H = H p ® T-L c 

Coin Evolution Operator: Any 2-dimensional uni- 
tary operator can be a coin evolution operator. The 
Hadamard operator is customary used in this work, 
where Uc is the Hadamard operator H: 



H 



1 



1 1 

1 -1 



1 



V2 



(|0)(0| + |0)(1| + |1)(0|-|1)(1|) 

Our system dynamics can be described by two steps, 
first, the coin operation, given by Uc €SU(2) acting only 
on 1-l c and the shift-position operation S which will move 
the particle accordingly to the state of the coin. There is 
a possibility of quantum superposition in these operators 
offering a phenomenon without a classical analogue. 

Conditional Shift Operator : S will move the par- 
ticle accordingly, transferring this way the quantum su- 
perposition to the total state in Ti. As with the previous 
operator, the only requirement is that of unitarity. A 
suitable conditional shift operator is given by: 

s = (j2 1* + wij ®i txti+fei*- 1>(^ ®u>u i, 

(12) 

The evolution of the system at each step of the walk can 
then be described by the total unitary operator: 



u = s{ip®u c ), 



(13) 



where Ip is the identity operator on T-L p and S and Uc are 
defined accordingly above. Note that if a measurement is 
performed after each step, we will revert to the classical 
random walk. 

Observables: Observables are defined in a standard 
way. In a standard way we simply mean the quantum 
mechanical form in which the action of an operator on 
its eigenstate gives an eigenvalue and they are what we 
call observables. We can either measure coin states first 
followed by position states or measure position states di- 
rectly. Results are the same because the operators com- 
mute. In the quantum walk we will of course not mea- 
sure the coin register during intermediate iterations, but 
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rather keep the quantum correlations between different 
positions and let them interfere in subsequent steps. The 
results of interference is what gives quantum walks its 
powerful behaviour over the classical case. 

Consider a particle with coin degree of freedom rep- 
resented in the spin. Consider a wave function rep- 
resentation of a spin-1/2 particle. The state is a 2-D 
vector |\&) = (1*0 )i |V^)) T > where the first part is the 
component of the wave-function of the particle in the 
spin- | f) state and the second one is the component 
in the spin-| \) state. The composite nature of the 
space is emphasized by using a tensor operator to ten- 
sor the space of the particle and the state is given by 



I*) 



= ™t 



1^ ) ® I t) + crHV^) ® I I); where we normalize 



the two wave-functions (i[>^\ip^) = (VHV^) = 1> so that 
1. The tensor product '®' separates the 



a 



t|2 



a 



l|2 



two degrees of freedom, spin and space, and will allow 
us to view the resulting correlations between these two 
degrees of freedom more clearly. The time development 
corresponding to a translation by I on the larger state- 
space of the spin-| particle can now be described by the 
unitary operator U. This operator induces a kind of con- 
ditional translation of the particle depending on its inter- 
nal spin-degree of freedom. In particular consider a situ- 
ation where the spin of the particle is initially in the state 
so that its wave-function is of the form \ip% ) ® | f)j 
then application of U transforms it to \i> Xo _j) ® | t) an d 
the particle will be shifted to the right by I. If the spin 
of the particle is in the state | |) , and so the total wave- 
function is given by \ip^ ) ® | t), then the translation 

operator will transform it to \ip^ +;) ® I t) an d the par- 
ticle will be shifted to the left. Quantum effects can be 
introduce and we have situations where the particle spin 
state is a cohenrent superposition of the two states. We 
can thus have an intial state of the form 



i 1 



IV%> ® (« T I t) 



(14) 



Application of the translation operator U will induce a 
superposition of positions 



U\*i 



- n t 



o^xo-l) ® 1 1) + a l \ip xa+l ) ® | 1). (15) 



In this work we have used a notations that are much 
related to the computational states. So in order to ensure 
a similar approach like in the simulation, we have used 
I = 1 and the position state of the particle is considered 
as |i) € Z. We would equally consider a situation only 
where the particle start in the porition i = and the 
coin state is state is a coherent superposition state of 

(|t> + U»A/2 



1. Hadamard walk 

Many forms of a coin operator can be used to but we 
consider discrete-time quantum walk is the walk using 
Hadamard operation as quantum coin operation and is 
known as the Hadamard wallP. A particle at origin in 



one of the basis state |0) or |1) of H c (internal state of 
the particle) is evolved into the superposition of H c with 
equal probability, by applying the Hadamard operation 



h = -L(|o)<o| + |o><i| 

such that 

(ip ® Uc){\ifo) ® 1 1» = 
tfp®#b)(hfo>®U» 



|1>(0|-|1><1|) 



(16) 



V2 
1 



KM ® (I t) + U» 



V2 



|Vo)®(|t)-U))-(17) 



The operation H is then followed by the conditional shift 
operation S in a general form, U = S(Ip ® Uc))- The 
evolution of the system is allowed to progress without 
intermediate measurements to evolve the particle in su- 
perposition of position space and realize a large number 
of steps of the Hadamard walk. After the first two it- 
erations of U , the left and right moving components of 
the amplitude begin to interfere, deviating from the clas- 
sical evolution, thus resulting in a quadratic speedup in 
the growth of the variance. The probability amplitude 
distribution arising from the iterated application of U is 
significantly different from the probability distribution of 
the classical random walk, as shown below. For the first 
two application, there is no difference between the classi- 
cal and the quantum case, but a difference is visible as the 
evolution gets to more than 3 steps. The particle initially 
in the state |0) drifts to the left and the particle with an 
initial state 1 1) drifts to the right as explained above using 
notion of state vectors. This asymmetry arises from the 
fact that the Hadamard operation treats the two states 
| f) an d | l) differently, phase difference of —1 in case of 
state | I). This phase difference, brings our some of the 
significant differences depending on the initial state of 
the particle contributes to the constructive interference 
on one side and to the destructive interference on the 
other side of the position space. Therefore, if the initial 
state of our particle is, for instance |0) (g> | f)i the first 
step of the quantum walk will be as follows: 



|0> ® | t) 



u. 



1 



v/2 

^ 2 
u 1 



= (|i)®|t) + |-i>®U». 

l(|2)®|t)-|0)®(|t)-4)) + |-2)®U)). 

(|3>®|t) + |l>®U> + |-l>®|t)-2|-l) 



2^2 



This illustrates that there is a similar probability to find 
the particle in position 0,1 and 2 for classical and quan- 
tum walks for zeroth, first and second steps. However 
there is a striking difference for the third and higher or- 
der, situation explained only by superposition principle, 
and this is the effect that is explored in most quantum 
application. Applying the coin unitary operator without 
intermediate measurements we can then obtain a plots 
shown below for below. The behaviour depends on the 
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nature of the initial state. When the coin initial state is 
not a coherent superposition, we obtain an asymmetric 
situation given by: The plot in figure [T] shows a plot for 



Probability 




Particle Position 



FIG. 1. Probability distribution of classical random walk 
(Red line Gaussian) and its quantum counterpart, quantum 
walk after N = 100 stepspl. The initial state is |0) ® | f). 



2. Discrete-time quantum walk on a periodic line 

The quantum walk on a particle in a periodic line can 
be simulated using an n-cycle. This can be viewed as 
a particle moving in a circle with nodes. This walk will 
then depend on the number or nodes in a circle or in other 
words the periodicity of the circle. Before introducing 
the concept of quantum limiting distribution, we provide 
an example of a quantum walk on a periodic graph: a 
discrete quantum walk on a cycle. A cycle is a graph G 
with vertex set = n|) and edges set E , G = (V, E). 

A quantum walk on G acts on a total Hilbcrt space H p <£> 
He . In this case of the Hadamard coin operator given 
above the shift operator is defined by S\0,j) = \0,j + 1 
mod n) and S\l,j) = |0,j — 1 mod n). The behaviour 
of quantum walk on a cycle as will be shown depends 
so strongly on the the number of nodes or periodicity of 
the cycle. Here we shall be concern with a circle only 
which is a cycle with \V\ = 2. Below is a a diagram 
showing quantum walk for a single particle on a circle 
with different number of edges or period T. 




-1 00 -80 



-40 -20 20 40 60 80 100 
Particle position 



FIG. 2. Probability distribution of classical random walk 
(Red line Gaussian) and quantum walk after each after N = 
100 stepsPl.The initial state is f |0) ® (| 1") + | 4). 




particle position 



FIG. 3. Quantum walk on different periodic cycle for a single 
walker. The probability distribution depends so much on the 
number of nodes and its shown for different number of nodes 
T. 



a quantum and a classical case after, N = 100 steps. The 
difference in the variance can been seen in the distribu- 
tion. The quantum case is shown for the non symmetric 
state |0) ® | f) • The asymmetry can be uplifted by by 
choosing the coin in a coherent superposition of the — f) 
and | 1} as shown below. 

The plot in figure [2] shows a plot for the classical and 
quantum walk after N = 100 time steps. Here the curve 
is symmetric because the coin start in a coherence super- 
position state. 



II. QUANTUM WALK OF TWO PARTICLE 

A. Introduction 

Quantum walk has been shown to have features with 
no classical analogueP^ . Just like in the previous sec- 
tion we will consider quantum walk in system with two 
particles. We consider two unitary operators to repre- 
sent our quantum equivalence of a coin and then we have 
two walkers. There are many features that can be ex- 
plored including the notion of superposition and entan- 
gled states^ for which no known classical equivalence 
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exist. Entanglement can be studied in the coin degree 
of freedom as well as in the particle degree of freedom. 
In this section we will present studies made only in the 
coin degree of freedom which we use an analogy of two 
particles with spin states that are entangled. Again we 
will present results that suggests that the behaviour of 
these walkers is also a function of the type of graph they 
traverse as we will compare and contrast evolution on a 
graph with and without periodic boundary conditions. 
In this section we will start by explaining quantum walk 
of two particle in a line and then we extend it to a cycle 
and make some comparisons to the results. 



B. Dynamics of two particle walk 

Two particles in a system can be interacting or non 
interacting. In the case of quantum walk, such a dynamic 
will live in a composite Hilbert space denoted by 



H12 =%i®U2 = H p ,i (g> H c ,i <8> H 



lp,2 



n 



c.2 



(19) 



where Hi and H2 represent the Hilbert spaces of parti- 
cles 1 and 2 respectively. Since the relevant degrees of 
freedom in our problem are the same for both particles, 
we have that both H\ and H.2 are isomorphic to H de- 
fined earlier for the one-particle case. The dynamics of 
of such a system is govern by a unitary given by 



U 12 = U ® U, 



(20) 



where U is given by equation ( 13 1 and is the same for 



both particles. Note also that in the case of identi- 
cal particles we have to restrict H12 to its symmetrical 
and antisymmetric subspaces, respectively for bosons and 
fermions. The other entangled states of the triplet states 
are not considered as all the of them are symmetric and 
gives the same results when considered. Let us then con- 
sider without lost of generality the case where both par- 
ticles start the quantum walk in the same position state 
|0,0), but with different coin states | \.) and | t)- Con- 
sider a simple way to understand this mathematically. 
Let us consider the coin subspace by using a representa- 
tion with |0) = I f) and |1) = | L). The coins consists of 
two qubits with basis {\j, k)},j, k € {0, 1}. The particle 
position space is given by {|m, n)}, ra, n € Z ,and then 
the generic state of the quantum walker at any time step 
t would be 



(21) 



We can also represent the evolution operator U = S(Ip <g> 
Uc), in this notation where, 

j,k j',k' 

is the coin operator, Ip is the identity matrix, and S 



is the shift operator given by equation 12 and which we 



represent here by equation [23] 

S\j,k)\m,n) = \3,k)\m + (-iy,n+(-l) k ). (23) 

We can see that the two particles both moves to the right 
when the coin is |0,0), and to the left when the coin is 
|1, 1). In the state |1,0), particle one moves to the left 
and particle two moves to the right, and the situation is 
reverse for the state |0, 1). 

Applying the evolution operator on state [21] we get 

Aj,k;m,n(t + 1) = ^ Cj,k;j' ,fc' Af M' !"»+(- \)i ,n+(- l) fe (*) ■ 
j>,k> 

,( 24 ) 

The probability distribution for the walker at position 
|ra, n) at time t is 



\ A ^;r. 

3,k 



(25) 



We can then go back to our normal formalism to rep- 
resent the possible states as is done below. There will be 
three cases we shall consider here. The pure separable 
state where our system's initial state will be given by: 



|(U)i|o,t) 2 . 



(26) 



Without any confusion, we have switched back to the 
normal notation where the '0' in the state above is show- 
ing the particle's position and has nothing to do with 
the coin degree of freedom. The other situation that is 
considered is for initial pure entangled states in the coin 
degree of freedom. The case of two maximally entan- 
gled coins is considered. The symmetric case where our 
system's initial state will be given by: 



1 

72~ v 



|^ +)i2 = ^(|0,|) 1 |0,t) 2 + |0,t)i|0,|> 2 ), (27) 



and the antisymmetric case where our system's initial 
state will be given by: 

|V>o~)i2 = i(|04)i|0,t)a - |0,t>i|<U>a). (28) 

These states differ only by the presence of a phase. These 
are states that are considered as bosonic and fermion par- 
ticles in physics when the notion of identical particles is 
introduced. After an evolution of N time steps a system 
starting with initial state |V ; o)i2 becomes; 



I^)i2 = ^|^')i2 = C/ iV |0,|) 1 C/ iV |0,t) 2 . (29) 



WVi 



N\ 



for initial conditions ( 26 ) . The probability distribution 
P 1 2(«, j; N) for finding particle 1 in position i and particle 



2 in position j for N = 100 steps is shown in section (III I 



We can show that the distributions of the particles are 
uncorrelated, and so P^ihj'i N) is simply the product of 
the two independent one-particle distributions: 



Pi2^r,N) 



P?(i;N)xPi(j;N), 



(30) 
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where Pf(i;N) is the probability distribution for find- 
ing particle 1 in position i after N steps, and similarly 
P^ij'jN) for particle 2. A corresponding results for 
movement on a cycle is given beside. 

In a similar fashion, when the particles are entangled, 
evolution of the system after N steps will be 



|V#>12 = #3 |V + )i2 = ^(0 N \Q,i) l U N \0.t) 2 +U N \OA),U N \Q:i) 2 



for the symmetric entangled state, and 



B. Absorbing boundary conditions. 

The following results illustrate the difference obtained 
for different initial states. The results shows some strik- 
ing difference. This is as a result of some nontrivial quan- 
tum correlations resulting from the measurements in the 
case of symmetric and antisymmetric states. The results 
is given below, 



\i>N)l2 = U N 2 \% 



> 12 



(u N \o,ihu"\OAh-u 



AT I 



V2 



for the antisymmetric case. 

The figures show the probability distribution for find- 
ing particle 1 in position i and particle 2 in position j in 
the "+" case, P&(i,j; N), and P^ihT, N) of the "-" case 
for N = 100. The effects of the entanglement are striking 
as noticed when comparing them as shown in the distri- 
butions plotted in section ( III I . There is a significant 



increase the probability of finding the particles in certain 
configurations and some are forbidden. This is one rea- 
son why there is a speed up in quantum algorithms based 
on quantum walk. The results for particle dynamics on 
a circle is also presented. There are some similarities as 
well as some striking differences. We can note from the 
graphs that the particles tend to finish as far as possible 
in the case of antisymmetric states a situation much com- 
mon to nature and also a situation that doesn't manifest 
itself well in periodic conditions, except for large enough 
periods compared to the total time time steps. To make 
a concrete comparison we will compare the dynamics of 
each particle for both the walk on a line and on a circle 
by considering their average positions and comparing for 
the three different states for both line and circle. We will 
also study the variations in their Von Neumann entropy 
for a general configuration and also for some specific po- 
sitions. 




FIG. 4. The probability distribution for two particles on a line 
after N = 100 time steps. The result is shown for product 
and symmetric states. 



Also in a similar way the results for product and anti- 
symmetric states is shown below. These behaviors are 




I | Position of panicle 1 



Position of panicle 2 



FIG. 5. The probability distribution for two particles on a line 
after JV = 100 time steps. The result is shown for product 
and symmetric states. 



III. COMPARISON OF RESULTS 

A. Introduction 

The discrete quantum walk of two particle on a straight 
line, shows interesting behaviours depending on the 
boundary condition. As shown irP^ there are many inter- 
esting behaviour that are a result of the initial state. The 
underlying dynamics of each particle is also presented, 
but in this chapter we shall compare results for vari- 
ous boundary conditions. We present results of quantum 
walk on a line for product, symmetric and antisymmetric 
coin state. 



investigated irP , in which they studied qualitatively 
the behavior of each of the particle using partial trace 
knowledge to get a reduced density operator for the evo- 
lution of each particle. Using our knowledge of bosons 
and fermions, we see that in the case of identical particles, 
in the symmetric state (bosons) , the particles tend to fin- 
ish together, while in the antisymmetric state (fermion), 
the particles tend to finish as far apart as possible, offer- 
ing a behaviour we expect, from their statistics. 



C. Periodic boundary conditions 

Using the results above a survey for the various ini- 
tial states, and then a implementation of the boundary 
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FIG. 6. The probability distribution for two particles on a line 
after N = 100 time steps. The result is shown for product, 
symmetric and antisymmetric states in absorbing and non 
absorbing boundary conditions. Here a Period T — 100 is 
used to ensure enough compatibility. 



conditions of a circle gives results that depends on the 
boundary conditions. A study made with two particles 
and for different types of boundary conditions is made. 
The probability distribution is periodic for smaller peri- 
ods. It can also be noticed from the distribution that 
some regions on the circle are favored. A comparison 
between periodic and absorbing boundary conditions is 
made below. There arc different ways to study the be- 
haviour of each of the situation. A reduced density could 
be used as already pointed above. One way to character- 
ize the behaviour of the three systems is to study their 
Von Neumann entropy as already pointed above. With 
the entropy of the combined system in mind, we can get 
a reduced description as well. With this we can study 
the behaviour of a single particle in the composite sys- 
tem. Various studies can then be made of the probability 
distribution for seeing the particle at position, i after k 
time steps. So fixing i the position we can obtain a plot 
for the probability distribution for seeing the particle at 
that position after each time step. 



D. Different periods 



A study of the period shows that the distribution de- 
pends so much on the periodicity of the system. Below 
are plots obtained for the symmetric state for the pe- 
riods T = 8, 16, 64 and 100. We make the plot in the 
form of a lattice, i.e a periodic lattice and the we can 
observe periodicity on the lattice as the particle undergo 
its dynamic. By lattice here, we mean that the particle 
after traversing the period, comes to a point that is iden- 
tical to the original pint and so we have a possibility to 
see the periodicity in the plot. Another way could be 
to make a plot that is limited to nodes equal to the pe- 
riod. This is less illustrative as it doesn't explicitly show 
the periodic nature of the results. A study of the prob- 
ability distribution to see a particle at position, i after 
k time steps, plotted for each i would can be made. A 
similar results for the one particle case is shown by ? in 
which the showed that the probability distribution for 
seeing the particle at the starting position is periodic up 
to T = 8 and that the system become really chaotic with 
an inability to predict where it can be from it entropy for 
T = 16, 32, 64 and 100. The results for different periods 



is shown figure III D 
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